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S: 

Abstract 

CN ■ We establish the soundness of the replica symmetric ansatz in- 

, troduced by M. Mezard and G. Parisi for minimum matching and 

I the traveling salesman problem in the pseudo-dimension d mean field 

0^ ' model for d > 1. The case d = 1 of minimum matching corresponds to 

. the 7r^/6-limit for the assignment problem established by D. Aldous in 

I 2001, and the analogous limit for the d = 1 case of TSP was recently 

■ obtained by the author with a different method. 

. We introduce a game-theoretical framework by which we prove the 

^ I correctness of the replica-cavity prediction of the corresponding limits 

also for d > 1. 



1 Introduction and background 



1.1 An example 

Suppose that the edges of the complete graph Kn for even n are assigned 
independent costs from uniform distribution on [0, 1], and that we study the 
minimum total cost of a perfect matching, in other words a set of n/2 edges 
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having no vertex in common. Based on the non-rigorous rephca method of 
statistical physics, M. Mezard and G. Parisi conjectured in [IB] that the min- 
imum cost converges in probabihty to '7r^/12. This was proved by D. Aldous 
[2], and it is known that the hmit depends only on the density of the cost- 
distribution at zero. Hence apart from a scaling factor, the 7r^/12-limit is 
stable under replacing uniform [0, 1] by various other distributions such as 
exponential, minimum of two independent [0, l]-variables etc. 

On the other hand there are distributions that are not of this form, but 
where the density at x > scales like a power of x as a; — > 0. Suppose 
for instance that the edge costs are sums of two independent uniform [0, 1]- 
variables. In that case [18] predicted that the minimum cost of a perfect 
matching is approximately 0.8086 ■ ^Jn. 

In this paper we obtain a rigorous confirmation of this prediction. We 
establish the correctness of the replica method for this and some related 
problems, notably the traveling salesman. We prove that the characteriza- 
tion in [18] of the limits (under rescaling by the appropriate power of n) is 
essentially correct. 

1.2 Background 

It has been known for some decades that methods of the statistical me- 
chanics of disordered systems apply to certain problems of combinatorial 
optimization. Much of the work in this direction stems directly or indirectly 
from G. Parisi's solution p6] of the Sherrington-Kirkpatrick model [B2| of 
spin glasses, established rigorously by M. Talagrand [35j. In S. Kirk- 
patrick and G. Toulouse suggested the mean field traveling salesman prob- 
lem (TSP) as an archetypal optimization problem sharing important features 
with spin glasses. M. Mezard and G. Parisi [ISl [ig EOl ED EZ] and Mezard 
and W. Krauth [T5] obtained several remarkably detailed predictions about 
minimum matching and the TSP with the replica and cavity methods. These 
predictions were based on the assumption of replica symmetry which is known 
to fail at low temperatures for models of spin glasses. It became clear that 
minimum matching and the TSP are different in this respect from models like 
the Sherrington-Kirkpatrick model and random /c-SAT. Several authors have 
verified the consistency of the replica symmetric ansatz by testing its various 
predictions numerically and theoretically [6l|7l[8lll0l[H[221l25l|28l[29l[30l|3l]. 

Replica symmetry is interesting from an algorithmic point of view since 
it is linked to the efficiency of Belief Propagation heuristics, see for instance 
[25] . Recently J. Salez and D. Shah [31J have obtained rigorous results in 
this direction for the assignment problem, and in view of our results, their 
conclusions should be valid also for the TSP. 
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In this paper we introduce a two-person game which is played on a graph 
with lengths associated to the edges. We show that on certain infinite graphs 
with random edge-lengths, this game has an almost surely well-defined game- 
theoretical value, and we argue that this property is essentially equivalent to 
replica symmetry for the minimum matching problem. 

The paper is organized as follows: In Section 12.11 we give further back- 
ground and state our main result on minimum matching. Theorem 12.11 The 
rest of Section [2] is devoted to the proof of this theorem. The most important 
new results are the introduction of the Graph Exploration game in Section [22] 
and the analysis of this game on the rf-PWIT in Sections I2.5H2.6I In Sec- 
tion [3] we state and prove the analogous theorem for the TSP, Theorem 13. 1[ 
The TSP corresponds to a comply- constrain version of Graph Exploration. 
We discuss some important differences in the analysis of this game compared 
to its normal counterpart, but avoid repeating arguments that are similar 
to those of Section [2j In Section H] we study the minimum edge cover. This 
problem has not been considered in the physics literature, but has a rich 
structure and can be analyzed by the same methods as matching and the 
TSP. Section O contains some concluding remarks. 

2 Minimum matching 

2.1 The mean field model and the minimum matching 
problem 

The mean field model of distance is a complete graph on n vertices whose 
edges are assigned i.i.d lengths (or costs) from a distribution on the pos- 
itive real numbers. If the edge-lengths are intended to model distances be- 
tween random points in d-dimensional space, then we expect P{1 < r) ~ c-r'^ 
for small r, since the probability of a point being within distance r of another 
should be proportional to the volume of a ball of radius r. Certain asymp- 
totical properties of optimization problems are known to depend (apart from 
trivial scaling factors) only on the parameter d, so that for instance uniform 
distribution on [0, 1] is equivalent to exponential distribution of mean 1, both 
belonging to the case d = 1. Similarly the sum and the maximum of two 
independent uniform [0, 1] variables both represent d = 2, and are equivalent 
apart from a scaling factor. 

Such results can often be established with standard techniques, but for 
convenience we make a specific choice of distribution for < d < oohj taking 
/ to be the d-th root of an exponential variable. To simplify the scaling in 
terms of n, we let / = (nXy^'^, where X is exponential of mean 1. In other 
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words / is the d-th. root of an exponential variable of mean n. This gives 

P{1 <r)= P{{nX)^^'^ <r) = P{X < r'^/n) = 1 - exp(-r'^/n) ~ r'^/n. 

We can thereby regard the lengths as generated from an underlying set 
of independent mean 1 exponential variables Xij. 

A favorite problem is minimum matching, which seems to be the simplest 
problem that allows the ideas of [TB] to be displayed in a nontrivial way. 
We ask for a set of edges of minimum total length such that each vertex is 
incident to exactly one. This obviously requires n to be even unless we allow 
one vertex to be left out of the pairing, but this is a minor issue since we are 
mainly interested in the large n asymptotics. 

The quantity of main interest is the total length M„ of the minimum 
matching. It is not difficult to guess roughly how M„ scales with n. From 
an arbitrary vertex, the order of the distance to the nearest neighbors is 
obtained by setting P{1 < r) ^ which leads to r ^ 1. If we believe 
that edge-lengths of order 1 will dominate the solution, then since a solution 
contains n/2 edges, we expect M„ to scale like n. 

It is natural to conjecture that Mn/n, which can be interpreted as the 
average cost per vertex of obtaining a matching, converges in probability to 
a constant depending on d. Our main result is that this is true for d > 1. 

Theorem 2.1. For every d> 1 there is a number ^M^d) such that 

^^hM). (1) 

n 

We believe that in principle the method applies also when < (i < 1, 
but we have run into some difficulties that have so far prevented us from 
establishing ([1]) in that case. 

Let us immediately state the easiest available bounds on /3m (c?)- For a 
lower bound we observe that the expected length of an edge in a matching 
must be at least the expectation of the length /min of the shortest edge from 
a given vertex. A factor 1/2 comes from averaging the total length of the 
n/2 edges over n vertices: 

Md) > 77^(Un) = - / P{lmin > v) dr = - exp{-r'' / n)^-' dr 

exp(-r'^)dr = -r(l + l/rf). (2) 

Getting an upper bound is not trivial, but it is known [2J that /3m (1) = vr^/12. 
For d> 1 the concavity of the mapping X 1— > (nXy^'^ gives the bound 

„ 1 /'"'^ 
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A calculation in fS] backed up by results in f24j gives the sharper bound 



1 p-Xf -X _ 1 _|_ 1 

valid for (i > (see also [35] and [SB])- A discussion of how to prove would 
be out of place here, but the idea is to use the matching that minimizes the 
sum of the underlying exponential variables Xij. 

It was observed in [37] that a greedy algorithm gives a matching of the 
right order of magnitude if > 1 (actually [2Z] considered the TSP, but 
leaving out every other edge of a tour obviously gives a matching). The 
greedy matching (also known as the Gale-Shapley matching) gives the weaker 
bound 

PM{d) < - ■ . / 

2 sm{TT/a) 

but is interesting in itself. 

Within the framework of the replica method, Mezard and Parisi obtained 
an analytical characterization of /5m (c?) which is conjectured to be correct for 
all > 0. They arrived at an integral equation which is equivalent to 

F{x) = exp (^-d l'^-^F{l - x) dl^ , (4) 

from which ^M^d) is obtained as 

PM{d) = ^ J J ix + y)'-'F{x)F{y)dxdy. 

—oo<x,y<oo 
x+y>0 

The method is inherently non-rigorous, and it has not been established that 
(jlj) has a unique solution except in the case d = 1, where the solution 

F{x) 



1 + 



leads to ^a/(1) = n^/l2. 

Mezard and Parisi also calculated /3m (2) numerically. In [IB], the edge- 
lengths have density rather than dl'^~^/n, which means that for d = 2, 
our edges are longer than those of [I8j by a factor ^J~rif2. On the other hand 
Mezard and Parisi considered a graph on 2N vertices but rescaled the total 
length of the matching by a power of A^, in this case ^JN rather than \/2N . 
Therefore the ground state energy E ^ 1.144 given in equation (24) of 
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is 2/?m(2). It is interesting to compare their numerical value to the bounds 
and ©. For d = 2 we get 

= 0.443 ■ ■ ■ < /5m(2) < • C(3/2) = 0.578 . . . , 

while |IS] gives Pa.i{2) ^ 0.572. The fact that the estimated true value is 
quite close to the upper bound indicates that the matching that minimizes 
the sum of the underlying exponential variables is not too far from the actual 
optimum. 

It is worth pointing out how the scaling works in a couple of simple 
examples. Suppose we take the distribution of edge lengths as / = max(f/, V), 
where U and V are independent and uniform in [0, 1]. Then 

P(l <r) = P{U < r)2 = 

if r < 1, and therefore the distribution belongs to the case d = 2. We 
determine the value of tq for which the expected number of points within 
distance tq from a given point is equal to 1. Approximately this happens 
when P{1 < To) ~ l/n, which gives 
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We can think of /3m (2) as the average cost per vertex of the minimum match- 
ing, measured with tq as the unit of length. Hence the total length of the 
minimum matching is approximately 

/3m (2) v^. 

If on the other hand we take the edge lengths to be distributed like U + V 
as in the introduction, then P{1 < r) ~ r^/2, and the unit of length is given 
by rQ/2 = 1/n or equivalently 

V2 

In that case the asymptotical total length of the minimum matching is 

Pm{2)V^. 

Although it does not follow from Theorem 12.11 our results apply also to 
the assignment problem, in other words minimum matching on the complete 
bipartite graph Kn,n- In the bipartite model the nearest neighbor distances 
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are the same and the only difference is that a matching contains n edges 
instead of n/2. With the two distributions max(f/, V) and U + V, the mini- 
mum assignments will have lengths approximately 2/5Af (2)y^ and /?M(2)v^8n 
respectively. 

A more precise value of /3m(2) was obtained by J. Houdayer, J. H. Boutet 
de Monvel and O. C. Martin [T3]. Using length-distributions normalized by 
the volume of the d-dimensional ball, they obtained the value 0.322580 for 
the limit. After sorting out the normalization, one finds that our /3m (2) is ^/^T 
times their value, which gives /?m(2) ~ 0.571758. Through an approximate 
solution of (jl]) we have obtained the value 

Pm{2) ^ 0.57175904959888. 

We have no theoretical estimate of the error, but these decimals seem to 
be stable. The values equivalent to /^^/((i) for integers d < 10 are given in 
Table 2 of [2] . Although we still do not know whether (j3]) has a solution, the 
numerical result can be regarded as an approximation of the fixed point of 
Ve (see Section 12.71) for an appropriately chosen 6, and therefore apart from 
the numerical error the result is backed up rigorously. 

On the mathematical side there has been considerable progress on the 
case d = 1. In particular the 7r^/6-limit in the assignment problem has 
received several different proofs [21 [IZl [211 EH] • From our point of view the 
result corresponds to the statement that for d = 1, the limit in ([1]) exists 
and /9m (1) = 7r^/12, but the asymptotic equivalence between assignment 
and matching on the complete graph is by no means trivial, and does not 
follow from [T71 [211 El]. The proofs together provide a quite detailed picture 
of the distribution of the total length as well as the local statistics of the 
optimum solution, and the analogous result for the TSP was established 
in |l0] . However, the proofs in pTl |2U [391 [10] are very different from the 
approach in the physics literature, and do not seem to generalize to d 1. 
The original proof by David Aldous [2] is the one that comes closest to 
justifying the replica symmetric ansatz (particularly in view of additional 
results in [31 [31]), but it seems to rely on finding a solution to ([1]). 

In the present paper, our aim is to show that the calculations in [18] are 
sound for quite general reasons. We prove that for > 1, ([T]) holds, and we 
characterize /3m (c^) analytically in terms of certain integral equations similar 
to dlj). Although we cannot find explicit solutions to these equations when 
d 1, our results show that the numerical computation of /?m(2) in [TS] is 
correct in principle. 

For d = 1, much more detailed results can be obtained. A more precise 
analysis of the d = 1 case with the present method, and a clarification of 
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its relation to the results of will be given in joint work with G. Parisi 
(manuscript in preparation). 

Our approach is "zero temperature", but similar to the replica-cavity 
method in that we reach the optimum solution through a limiting process. 
We introduce a parameter and study "diluted" problems where partial 
matchings are allowed but penalized hy 9/2 for each unmatched vertex. The 
original problem is recovered in the limit 9 ^ oo. The parameter 9 plays a 
role similar to the inverse temperature in statistical physics. Finite 9 allows 
for a certain local freedom that destroys all long-range interactions. In par- 
ticular, adding or deleting a vertex has only a local effect on the optimum 
solution. In [18] a similar assumption seems to be crucial for the renormal- 
ization that leads to (jl]). 

2.2 Graph Exploration 

The following two-person zero-sum game was invented in an attempt to find 
a mathematically sound interpretation of (jlj). We call it Graph Exploration 
since it somehow centers around the question whether it is worth the price 
to be the first to explore a new part of the graph. We are given a graph with 
nonnegative edge lengths, a starting point and a nonnegative parameter 
9. Alice and Bob take turns choosing the next edge of a self-avoiding walk, 
with Alice starting the game from v. The player who makes a move pays the 
length of the edge to the opponent. At each turn, the moving player also has 
the option to, instead of moving, terminate the game by paying 9/2 to the 
opponent. Each player tries to maximize their total payoff. 

Notice that there is no randomness in the game. The players are assumed 
to have perfect information about the graph including the edge-lengths. We 
can immediately make some observations: 

• If the graph is finite, then there is a well-defined game-theoretical value. 

• If the graph is infinite, there may or may not be such a value. For 
instance, if all edges have the same length / < 9^ then no player will 
ever want to terminate the game. 

• Edges of length more than 9 are irrelevant to the game. If Alice moves 
along such an edge, then Bob can terminate the game, and even though 
this may not be Bob's best option, it would still have been better for 
Alice to terminate in the first place. 
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2.3 The diluted matching problem 

There is a relaxation of the minimum matching problem that we refer to as 
the diluted matching problem. Instead of requiring each vertex to be covered 
by the matching, we allow for any partial matching, with a penalty of 6/2 
for each vertex that is not matched. 

For the moment we regard the parameter 6 as fixed. If G is a finite graph 
with given edge lengths, we let M{G) be the cost of the diluted match- 
ing problem. More precisely, M{G) is the minimum, taken over all partial 
matchings, of the sum of the edge lengths in the matching plus 6/2 times the 
number of unmatched vertices. 

Proposition 2.2. Let G be a finite graph with given edge lengths, and let v 
be a vertex of G chosen as the starting point for Graph Exploration. Then 
Bob 's payoff under optimal play is 

M{G)-M{G-v). 

Proof. Suppose that the neighbors of v are f i, . . . .Vk, and that the edges from 
V to these neighbors have lengths li, . . . Jk- Let f{G, v) be Bob's payoff when 
the game starts at the vertex v. By minimizing over Alice's move options, 
we recursively characterize / by 

fiG,v) = mmi6/2, k - f(G - v,v,)). 

On the other hand, the cost of the diluted matching problem satisfies 

M{G) = mm{6/2 + M{G -v),k + M{G -v- Vi)). 

Subtracting M(G — v) from both sides, we see that 

M{G) - M{G - f ) = min(^/2, k - {M{G - v) - M{G -v- Vi))). 

This shows that f{G,v) and M{G) — M{G — v) satisfy the same recursion, 
and it follows by induction that they are equal. □ 

It is clear from Proposition 12.21 and its proof that Alice achieves opti- 
mal payoff by starting along the edge of the optimal diluted matching, if 
there is such an edge from v, and otherwise by terminating immediately. By 
induction it follows that consistently playing along edges of the optimum di- 
luted matching, and terminating when no such edge is available, is minimax 
optimal. Therefore under mutual optimal play, the path described by the 
game is the symmetric difference of the optimal diluted matchings on G and 
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G — V. Actually the argument provides a simple proof of the fact that this 
symmetric difference is a path. 

Since the diluted matching problem can be solved efficiently by standard 
matching algorithms, it follows that Graph Exploration can be played op- 
timally with a polynomial time algorithm, but from our perspective this is 
beside the point. The advantage of introducing the game is that if the graph 
is infinite there may still be a well-defined game-theoretical value. This value 
then replaces M{G) — M[G — v) and allows for the equivalent of the renor- 
malization argument of [18] in a mathematically consistent way. 

2.4 Approximation by the PWIT 

The Poisson Weighted Infinite Tree (PWIT) was introduced by Aldous [H [2] . 
The PWIT is a rooted tree where each vertex has a countably infinite set 
of children, and the edges to these children are assigned lengths given by 
a rate 1 Poisson point process on the positive real numbers (independent 
processes for all vertices). The PWIT is a local weak limit of the mean field 
model, a statement which has been made precise in slightly different ways 
in the literature. We establish a simple version which is convenient for our 
purpose. We first treat the case d = 1, and later establish an easy refinement 
valid for general d. Recall that for d = 1, the edges in the mean field model 
are exponential of mean n. 

For given 6 and a positive integer k, the {k, 6) -neighborhood of a vertex 
f in a graph is the subgraph defined as the union of all paths from v of at 
most k edges, each of length at most 9. 

Lemma 2.3. Suppose k and 9 are given. Consider the graph Kn with a 
specified root v, and a random process consisting in assigning independent 
exponential lengths of mean n to the edges. There is a coupling of this process 
to the PWIT such that with probability at least 

(g + 2)fe 

^1/3 ' 

the [k, 9) -neighborhoods of v and of the root of the PWIT are isomorphic, 
with corresponding edges having equal length. 

Proof. We start from a PWIT rooted in a vertex v', and assign lengths to 
the edges of Kn through a random mapping of the {k, ^)-neighborhood of v' 
to Kn- We start by mapping v' to v. Then we sequentially map the vertices 
of the {k, ^)-neighborhood of v' to independent uniformly chosen vertices of 
Kn through a tree search (say depth-first). If we ever choose the root f or a 
vertex that has already been chosen, then we let the procedure fail. 
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To see that this is compatible with the probabihty measure on the edge 
lengths of Kn-, define an extended model in the following way: For each pair 
of vertices in Kn there is an infinite sequence of edges whose lengths are 
given by a rate 1/n Poisson process on the positive reals. Moreover, for each 
vertex there is a sequence of loops whose lengths are also given by a rate 1/n 
process (hence this model differs slightly from the friendly model of pO]). 
The original model is then recovered by discarding all loops and all edges 
except the shortest one between each pair of vertices. Now we explore the 
(/c, ^)-neighborhood of v in the extended model through a depth-first search 
from V, and "fail" if that neighborhood turns out not to be a tree. 

We want to estimate the probability of failure. Let be the number of 
vertices in the (/c, ^)-neighborhood of v' (including v'). Then 

E{N) = i + e + e'^ + --- + e^. 

Conditioning on A^, the expected number of collisions is 

n n 

Therefore the probability of at least one collision is at most N^/n. Now 
P(failure) < P(failure | < n^^^) + P{N > n^^^) 

^ ^ 1 + e + + ■■■ + 6'' _ 2 + e + e^ + --- + e^ ^ (g + 2)^ 

— ^ ^1/3 ^1/3 — j^l/3 

□ 

Lemma [23] can easily be generalized to the neighborhoods of several ver- 
tices, with the same method of proof: 

Lemma 2.4. Suppose m vertices in Kn are chosen independently of the edge 
lengths. Then with probability at least 

{mO + m + 1)^ 

the union of their {k^ 9) -neighborhoods is isomorphic to a disjoint union of 
the {k, 9) -neighborhoods of the roots of m independent PWITs. 

For general d, we get the edge lengths by raising the underlying exponen- 
tial variables to the power 1/d. To obtain a coupling, we introduce the 
d-PWIT, which is just the ordinary PWIT modified by raising the edge 
lengths to power 1/d. The original PWIT is the 1-PWIT, and the {k,9)- 
neighborhood of the root of the 1-PWIT corresponds to the (k, ^^/'^)-neighbor- 
hood of the root of the d-PWIT. 

By rescaling, the generalization of Lemma 12.41 becomes: 
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Lemma 2.5. Let k, 6 > 0, n and d > 1 be given. Consider the pseudo- 
dimension d mean field model on n vertices, with m vertices Vi, . . . ,Vm chosen 
independently of the edge-lengths. 

There is a coupling of this process to m independent d-PWITs such that 
with probability at least 

{mO'^ -\- m + 1)'' 

the union of the {k, 6) -neighborhoods ofvi,...,Vm is isomorphic to the {k,6)- 
neighborhoods of the roots of the d-PWITs, with corresponding edges having 
equal length. 

2.5 Graph Exploration on the d-PWIT 

Here we assume that d > 1, although some of the results hold also for 
< d < 1. In view of the results of Section 12.41 it makes sense to study 
Graph Exploration played on the d-PWlT. If f is a vertex of the d-PWlT 
we let TqIv) be the subgraph that can be reached from v by downward paths 
consisting of edges of length at most 9. The subgraph Tg{root) is called the 
6-cluster, and clearly nothing outside the ^-cluster is relevant for the game. 
Notice that the underlying graph of the ^-cluster is a Galton- Watson tree 
with Poisson(6''^)-distributed offspring. 

Our main objective is to show that although a priori the game does not 
need to terminate, there is almost surely a unique sensible way of assigning to 
it a game-theoretical value. The precise statement is Proposition 12. 101 below, 
and this is the key to the proof of Theorem 12.11 Proposition 12.101 shows 
that as ^ oo, a certain form of symmetry-breaking does not occur. This 
appears to be the fundamental reason why the replica and cavity methods 
are correct for the matching problem. 

When we speak of the value of a vertex v, by convention we mean the 
value of having moved to f , in other words the value of playing second if 
the game was played on TqIv) starting from v. If such a value f{v) can be 
defined consistently, it must clearly satisfy 

f{v)=mmi9/2,k-f{u:)), (6) 

where k is the length of the edge to the r.th child Vi of v, and the minimum 
is taken over 6/2 and the sequence of /j — /(f «) as Vi ranges over all children. 

For a given realization of the ^-cluster, we say that a function / from its 
vertices to the real numbers is a valuation if it satisfies A valuation can 
be regarded as a consistent way for a player to assess the positions of the 
game. We observe the following: 
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• A valuation must satisfy — 6'/2 < f{v) < 6/2 for every v. 

• If f is a leaf of the ^-cluster, then f{v) = 6/2. 

• If the ^-cluster is finite, there is a unique valuation. 

Proposition 2.6. There is almost surely a valuation. 

The only reason we say "almost surely" is that we haven't excluded the 
possibility that a vertex may have infinitely many children in the 6'-cluster. 
If this was the case, we would have to replace minimum by infimum in 
but this is an event of zero probability. 

Proof. Consider a "partial valuation" obtained by assigning values in 
favor of Bob to the vertices at distance k from the root. More precisely, 
these vertices get value ^/2 if is even and —6/2 if k is odd. Values are then 
propagated towards the root according to ([6]). As k increases, the values 
/^(f ) form a monotone sequence at each vertex v (decreasing at even levels, 
increasing at odd levels). Therefore there is a pointwise limit 

feiv) = lim f%{v), 

fe— >oo 

and it is easily verified that Jb is a valuation. □ 

Clearly Jb is at least as favorable to Bob as any other valuation. We 
can order the valuations from Bob's point of view by saying that /i < /2 if 
whenever v is at even distance from the root, /i(f ) < f2{v), and whenever v 
is at odd distance from the root, /i(f) > f2{v). Under this ordering the set 
of valuations forms a lattice where Jb is the maximal element, and similarly 
there is a minimal element which is most favorable from Alice's point of 
view. 

We are aiming to show that almost surely Ja = fs- This holds trivially 
in the range 6 < 1, since the 6'-cluster is almost surely finite. For ^ > 1, 
the 6'-cluster is infinite with positive probability, and the scenario that we 
wish to exclude is that at some critical value of 6 there occurs a breaking of 
symmetry after which J'a is distinct from fs- 

The question whether = fs is in a curious way similar to questions 
of the efficiency of game-tree search in games of perfect information such as 
chess. Uniqueness of valuation means that a game-tree search will be effec- 
tive, while symmetry-breaking corresponds to a situation where important 
long-term features of a position stay invisible to any fixed-depth search. 
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2.6 The branching of near-optimal play 

For the moment we take fs as our default valuation. This defines a strategy 
in an obvious way: From a vertex v, terminate if /^(f ) =6/2, and otherwise 
move to the child Vi for which fsiv) = U — f{vi). There seems to be the 
possibility of a tie in which several move options would be consistent with 
/b, but /b has the property that fsi'Vi) depends only on TQ{vi). Therefore 
h — fsivi) has continuous distribution and is independent of Ij — fB^Vj) for 
i j- It follows that the probability of a tie between move options is zero. 

Let 6 > 0. We say that a move from v to Vi is optimal if U—fBivi) = /^(f ), 
and 6-reasonable if Zj — fsivi) < fsiv) + 6. Let R be the subtree of the 9- 
cluster formed by all paths from the root consisting of 5-reasonable moves by 
Alice and optimal moves by Bob (a move can be 5-reasonable even if Zj > ^ 
so some 5-reasonable moves are excluded, but this is not important). Let 
R{k) be the set of vertices of R at distance k from the root. 

Proposition 2.7. If 6 is sufficiently small, then R is almost surely finite. 

We let 

H{k) = # {i; e R{k) : fB{v) < 6/2} + \-i^{ve R{k) : fB{v) = 6/2} . (7) 

The proof of Proposition 12.71 consists in showing that for sufficiently small S, 
EH{k) -> as /c ^ oo. 

The event v G R{k) does not depend on TqIv) through anything else than 
fB{v). It follows that if we condition on t> G R{k) and on fsiv), the structure 
of Tq{v) is distributed as if we condition on fB{v) only. Therefore we first 
assume that v is an "arbitrary" vertex of the 6'-cluster in the sense that Tg{v) 
is itself equal in distribution to Telroot). The children of v in the 6'-cluster 
are denoted by Vi. 

Lemma 2.8. The points (Zj, fsivi)) constitute a two-dimensional inhomoge- 
neous Poisson point process on the square [0,^^] x [—6/2,6/2]. 

Proof. The sequence of edge lengths U is a Poisson point process. Since 
fB^Vi) depends only on Te^Vi), the /B(fi)'s are independent of each other 
and of the Zj's. □ 

By the I- f -square we mean the square [0,6^] x [—6/2,6/2]. We let /x^ be 
the measure on the Z- /-square associated with the Poisson process of pairs 
(Zj, fsivi)). The measure is degenerate on the line / = 6/2 in the sense that 
this line has positive measure. Also notice that we do not assume that these 
point processes are equal in distribution for all v (this is what we are about 
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to prove). Prom what we have estabhshed so far it is conceivable that /ly 

depends on whether v is at even or odd distance from the root. 

To bound EH(k: + 1) in terms of EH{k) we bound the expected number 
of moves in R from a vertex v G R{k) in four cases, depending on whether 
Ahce or Bob is about to move and conditioning either on /^(f ) < 9/2 or on 
/s('f) — ^/2. The calculations rely crucially on the fact that for d > 1 the 
density of I is increasing, and that therefore the measure of a subset of 
the /-/-square increases under translation to the right. 

We first consider the case that Alice is about to move from a vertex 
V e -R(/c), where thus k is even. Suppose first that /^(f) < 0/2. Alice's 
optimal move is given by a point (Zj, fB{vi)) above the diagonal I — f — 9/2 
in the /-/-square. If we condition on fsiv) G [a,b] for some a,b such that 
-e/2 <a<b< 9/2, then 

^[jB(,Vi) - 9 2) - — 

livil - / e [a,b\) 

Hy{a + 9/2<l<b + 9/2) f^M - 9/2) > eM'O )■ (8) 

It follows that the probability that fnivi) = 9/2 conditioning on v G R(k) 
and fsiv) < 9/2 is at least exp(— ^^°'), and that therefore the optimal move 
by Alice contributes to EH{k + 1) by at most 1 — 1/2 • exp(— 6''^). 

The expected number of non-optimal 5-reasonable moves is at most 5 ■ 
d9^~^ = o(l) as 5 — > 0. Hence the expected contribution to if (A; -|- 1) when 
Alice moves from a vertex v such that fsiv) < 9/2 is at most 

l-lexp(-^'^)+o(l). 

By o(l) we mean a term that can be made as small as we please by making 
S small. 

Consider now the case that Alice moves from a vertex v e R{k) with 
fsiv) = 9/2. Then there is no optimal move (the optimal decision is to 
terminate), and again the expected number of 5-reasonable moves is at most 
5 ■ d9'^-^ = o(l). It follows that 

EH{k + 1) < ("i _ 1 exp(-^'^) + 0(1), <l-- + o(l). 



EH{k) - V 2 ' ^ " 1/2 J - 2 

When Bob moves, there is no optimal move if /b(w) =9/2 and at most one 
if fsiv) < 9/2. Hence the growth factor for H{k) over a pair of moves, one 
by Alice and one by Bob, satisfies 

EHik + 2) 1 , . 

^L_i<l--expM<') +0(1X1, 
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uniformly in if 5 is sufficiently small. It follows that EH{k) ^ as A; ^ oo 
and this completes the proof of Proposition 12.71 

The upper bounds on the expected contributions to H{k+1) when moving 
from a vertex v G R{k) are summarized in the following table: 



Player to move 




Vertex 


Contribution to 


EH{k + l) 


Alice moves 


fB 


(v) 


< e/2 


1 - 1/2 • exp(- 


-e'') + o{i) 




fB 


{v) 


= e/2 


o(l) 




Bob moves 


fB 


{v) 


< e/2 


1 






fB 


(v) 


= e/2 








Lemma 2.9. For sufficiently small 6, there is almost surely no infinite path 
starting anywhere in the e-cluster and consisting of optimal moves by Bob 
and 6-reasonable moves by Alice. 

Proof. If such a path started from the root, it would be a subset of R, and 
R is almost surely finite. This event therefore has probability zero, and it 
follows that the probability of such a path anywhere in the 6'-cluster is also 
zero. □ 

Proposition 2.10. There is almost surely only one valuation. 

Proof. It suffices to show that almost surely fA{root) = fB{root). Suppose 
therefore that this is not the case. Now let both Alice and Bob play "op- 
timistically" in the sense that Alice plays according to fA and Bob plays 
according to fs- Obviously they can never agree on an outcome of the game, 
so play has to continue forever. From Bob's perspective, it will seem that 
Alice sometimes makes mistakes that improve Bob's position. On the other 
hand the total gain (from Bob's perspective) of all these mistakes cannot be 
more than e, because the moment it adds up to more. Bob can terminate the 
game and receive a payoff greater than /^(root), and thereby also greater 
than /yi(root), which is a contradiction. Therefore the game must eventually 
reach a point where Alice's all future mistakes relative to add up to at 
most 6. The play from that point on will contradict Lemma 12. 9[ □ 

We need no longer distinguish between and fs, and we denote the 
almost surely unique valuation by /. Now recall the partial valuations 
and define similarly by choosing the values at level k in favor of Alice. 
Notice that and f\ are the upper and lower bounds on / that we get by 
looking k moves ahead from the root. 

Proposition 2.11. E (^f^{root) — /^(root)) — > as k ^ oo. 
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Proof. We have established that almost surely there is only one valuation. 
This means that almost surely, /^(root) — f\{root) monotonely as A; — > 
oo. The statement now follows from the principle of monotone convergence. 

□ 



2.7 Interpretation in terms of integral equations 

We want to obtain, to the extent possible, an analytical characterization of 
the distributions of /^(root), f^{root), and their common limit f{root). We 
have 

f\+\root) = min(^/2, k - /^'(^O), (9) 
where Vi ranges over the children of the root. Notice that 

Clearly the same holds with the roles of Ahce and Bob interchanged. 

Suppose now that we describe the distribution of f%{root) by the function 

Gk{x)=P{f%{Toot)>x), 

and similarly 

Fu+i{x) = P{f\^\root)>x). 

Then for —9/2 < x < Q jl^ is the probability that there is no event 

in the inhomogeneous Poisson process of Vi such that /j — f)^^{vi) < x, or 
equivalently, that there is no li such that f^^{vi) > li — x. Here it doesn't 
matter whether the inequality is strict or not, so for given x and li, 

P{ft\v,)>h-x)=Gkih-x). 

The sequence of li such that fA^^{vi) > k — x is therefore the set of points 
in a thinned Poisson point process of rate dl^~^Gk{l — x), and it follows that 

Fk+i{x) = exp ( -d J l'^-^Gk{l - x) dl 

Therefore we define an operator Vg on functions on the interval [—9/2,9/2] 

by 

r'e/2+x 

{VeF){x) = exp I -d I f-^F{l - x) dl 



We have F^+i = VeiGk), and by reversing the roles of Alice and Bob, Gk+i — 
Ve{Fk). The distributions of f\{root) and f^{root) are thus obtained by 
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starting from Fq = and Gq = 1 (on the interval [—6/2,6/2]) and iterating 
the operator Vg. But since Gi = Gq it follows inductively that F2 = Fi, 
G3 = G2 and so on. Therefore in reality there is only one sequence of 
functions, obtained by iterating Vg starting from the zero function. 

The operator Ve is decreasing in the sense that if F{x) < G{x) for every 
X, then {VeF){x) > {VgG){x) for every x. It follows that if we start from the 
function which is identically zero (or identically 1) and iterate, the sequence 
of functions must either converge to a fixed point or approach an attractor 
of period 2. Proposition 12.101 is equivalent to the statement that for every 
6 > and every d > 1, the sequence converges to a fixed point. Actually it 
is easy to see that if we start from any real integrable function F, then after 
two iterations we have a function which takes values in [0, 1], in other words 
lies between Fq and Gq. Therefore the subsequent iterates will be squeezed 
between Fk and Gk and thus converge to the same fixed point. In particular 
Vg has only one fixed point. 

The similarity to the Mezard-Parisi integral equation (jl]) is clearly visible. 
Naturally we may define an operator V^o by 



It seems clear both from numerical evidence and in view of the results 
we have established, that as 6^ ^ cxd, the fixed point of Vg should converge 
uniformly to a limit function which is a unique fixed point to Voo, in other 
words a unique solution to the Mezard-Parisi equation (jl]). We certainly 
believe that a more detailed analysis will show this to be true (possibly the 
ideas of [31] can be extended to c? > 1), but we leave it as an open conjecture 
since it is not necessary for our proof of Theorem 12. 1[ Moreover, the natural 
way to obtain numerical results from (jl]) is to approximate F{x) by 1 for 
large negative x and by for large positive x. Therefore in practice the 
numerical results based on (jl]) reduce to to solving the equation Vg{F) = F 
on a bounded interval. 

2.8 The density of the minimum diluted matching 

We now return to the mean field model Kn on n vertices. Suppose that 
6 and d > 1 are fixed and let the random variable g„ be the proportion 
of vertices that are not matched (for which we pay the punishment of 6/2) 
in the optimum diluted matching. Here and in the following we let q = 
P{f = 6/2) = F{6/2), where F is the fixed point of Vg. In other words q 
is the probability that Alice quits immediately in Graph Exploration on the 




d-PWIT. 



18 



Proposition 2.12. As n ^ oo, q. 

Proof. We show that Eqn — > q and var(g„) — > 0. Let A; be a positive integer. 
With probability 1 — o(l) as — > oo, the {k, 6')-neighborhood of a given 
vertex v in is isomorphic to the first k levels of a d-PWIT. 

By choosing large k, we can make E (^f^{root) — f\{root)) as small as we 
please, and provided that the coupling to the d-PWIT succeeds, the game the- 
oretical value of Graph Exploration on Kn starting at v is between f\{root) 
and f^iroot). Therefore conditioning on success of the coupling to the d- 
PWIT, 

P{f\{root) = e/2) < P{v is not matched) < P{fl{root) = 9/2), 

and both sides converge to g as A; — > oo. 

To bound the variance of g„ we simply take two vertices Vi and V2 of Kn 
and estimate the probability that neither is matched. To do this we apply 
Lemma 12.51 with m = 2. With high probability the {k, 6')-neighborhoods of 
f 1 and f 2 are disjoint and isomorphic to the {k, ^)-neighborhoods of the roots 
of two independent ci-PWITs. It follows that the probability that neither is 
matched converges to g^. □ 



2.9 The cost of the minimum diluted matching 

We wish to find the normalized limit cost of the minimum diluted matching. 
This cost splits naturally into the length of the participating edges and the 
cost of the penalties for the unmatched vertices. The penalties have been 
taken care of in the previous section, and therefore we concentrate on the 
participating edges. We let Mn{9) be the total length of the participating 
edges in the optimum diluted matching. 

Theorem 2.13. For each 9 and d>l, there is a number I3e{d) such that 

^^AW. (10) 

n 

Proof. Recall that the edge lengths are distributed like (nX)^/'^, where X 
is exponential of mean n. Therefore the density function of the length of a 
single edge is 

dl'^^^ ( 

exp . 

n \ n ) 
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The expectation of Mn{9) is the total number of edges in the graph times 
the expected contribution to M„(6') from a single edge: 

EMn{9) = Q 'n 'j(^ ^'^ ' ^^P(~^'^/'^) ■ ^(participation given length /) dl. 

(11) 

Deleting the factor exp{—l'^/n) will introduce an error of at most a factor 
(1 — 6'^/n). Normalizing to obtain a quantity of order 1, we get 



EM„{9) _ d 
n 2 Jo 



l'^ ■ P(participation given length I) dl + o{l). (12) 



We now choose a positive integer k. We explore the {k, 6')-neighborhood 
of the endpoints u and v of the edge e and discard the cases of "failure" when 
we cannot successfully couple to two independent d-PWITs to depth k. By 
choosing k suitably as a function of n we can make k tend to infinity while 
the probability of failure is o(l). 

Given that the coupling succeeds, the maximum length at which e par- 
ticipates lies between f^iu') + f\{y') and f^iu') + /^(f ')> where u' and v' are 
the roots of the two d-PWITs. Replacing these bounds by f{u') + f{v') will 
introduce another error of o(l). Hence f|T2l) is equal to 



\- j\''-P{l<f{u') + f{v')) dl + o{l). 

Here f{u') and f{v') are independent and satisfy P{f > x) = F[x) where F 
is the fixed point of Vq. By partial integration it follows that 

y- j j {x + yY-^F{x)F{y)dxdy. (13) 



EMn{e) 



-e/2<x,y<e/2 
a;+j/>0 



We denote the right hand side of (fT3l) by Pei^d). To see that (fT3!) can be 
strengthened to convergence in probability as stated in (ITU]) we again apply 
Lemma 12.51 this time with m = 4. It follows that the expected contribution 
from an arbitrary pair of edges to the square of Mn{9) is asymptotically 
the same as the square of the expected contribution of one edge, and that 
therefore var(M„(6')) = oin^). □ 



2.10 Perfect matching 

Here we complete the proof of Theorem 12.11 Assuming that n is even, we 
study the length M„ of the minimum perfect matching. Naturally we expect 
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perfect matching to correspond to infinite 6, and the remaining step essen- 
tially amounts to showing that we can interchange the order in which n and 
9 go to infinity. For the bipartite graph and d = 1 this was proved in [1]. 
Without claims of originality we give a self-contained proof valid for d > 0. 
This proof is based on expander properties of random graphs along the same 
lines as [TT]. I thank David Aldous for pointing out that the method of [TTj 
applies here. 

Recall that Pe{d) is defined as the right hand side of f|T3l) . Clearly l3e{d) 
is upper-bounded according to ([3]) and increasing in 9 by ffTOj) . We define 



as n — > oo. 

For every e > and every g > 0, we can find an n and a 9 such that with 
as high probability as we please, the minimum diluted matching has cost at 
most {PM^d) + e) ■ n and density at least 1 — g in terms of vertices covered. 
Therefore in order to complete the proof of Theorem 12.11 it suffices to show 
the following: 

Proposition 2.14. With high probability, a partial matching that covers a 
(1 — q) -fraction of the vertices can be completed to a perfect matching in a 
way that increases the total length by at most 6n, where 6 depends on q but 
not on n, and 6^0 as q —>■ 0. 

We introduce some extra edges by letting each pair of vertices give rise 
to a Poisson process of edges. More precisely, for each pair of vertices there 
is a rate 1 Poisson point process on the positive real numbers, and we let the 
sequence of edges have lengths given by (nXiY^^ where Xi are the points of 
the process. Obviously the extra edges do not change the minimum matching. 
We randomly color every edge red or green, where the probability of red is 
1 — p and the probability of green is p, for some p which will go to zero as n 
goes to infinity. 

First we find the minimum diluted matching M^ed on the red edges. Then 
before looking at the green edges we choose arbitrarily a bipartition of the 
vertices into two sets A and B of size m = n/2 such that every edge of M^ed 
connects a vertex of A to a vertex of B. Then we look at the green edges that 
connect A to B and give each of them a random orientation by independent 
coin flips. We let D be the set consisting of the 13 cheapest green edges 
directed from each vertex to the opposite side of the partition. 



Pnid) = lim Pe{d). 

9—fOo 

What remains is to show that for every e > 0, 
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Lemma 2.15. With high probability D has the following expander property. ■ 
If S is a set of vertices from one side of the partition, and 1 < l^l < m/3, 
then \S'\ > 2 15*1, where S' denotes the set of D -neighbors of S . 

Proof. If this condition is violated, then there is a positive integer s < m/3 
and a set of s vertices on one side of the partition such that all its 13s edges 
go into a certain set of 2s vertices on the other side. The probability that 
this happens is at most 

- E (:)(:) (I)'" 

l<s<m/3 \ / \ / \ / 

Using the standard inequality 




we find that the failure probability is at most 

2- V e^'2^^'(—Y^' 
^ \m) 

l<s<m/3 

By log-convexity of the summand, the maximum of the terms is attained 
by either the first or the last one. The first term (s = 1) is 0(l/ri^°) and the 
last one is at most 

/„3oll\ W3 

(^) 

For large n the first term will dominate, and since there are Oin) terms, the 
failure probability of the expander property is 0(l/n^). □ 

The number 13 was chosen to make the constant in the left hand side 
smaller than 1. If the expander property does not hold, then our scheme 
fails. If we want to keep the expected cost of completing the matching to 
oin) we may have to replace the number 13 by a larger number depending 
on d. We can then, in the cases of failure, just pick an arbitrary matching of 
red edges and the expected length will be 0{rt}^^l'^\ This will do provided 
that the probability of failure is sufficiently small. 

Proof of Proposition \2.14\ If the expander property holds, then we complete 
Mjied to a perfect matching by using the edges in D. We extend the matching 
successively by finding an alternating path that connects two unmatched 
vertices. When k unmatched vertices remain on each side, we can find such 
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an alternating path of length 0(log(n/fc)). A simple calculation shows that 
the total number of edges of D that become involved in completing the 
matching is 0{nq\og{n/ q)). Notice that this bound is deterministic and 
holds whenever the expander property holds. 

The completion of the matching is done independently of the actual 
lengths of the edges in Z^, so the expected total cost of completing the match- 
ing (given that the expander property holds) is bounded by 0{nq\og{n/ q)) 
times the expected length of an edge in D. Now it is clear that if p as 
— > cxD, we can keep the total length of Mrscl to Mn{6) + o{n) while at the 
same time keeping the cost of the completion process to o(n). □ 

This completes the proof of Theorem 12.11 



2.11 The 7rVl2-limit for d = 1 

The case d = 1 corresponds to the model studied by Aldous in [H, [2]. In 
our terminology his result (anticipated in [IB]) is that /3m(1) = 7r^/12. We 
briefly show how to derive this from our present approach. In [21 [18] the 
calculations start from equation (jl]) which has the solution F{x) = l/(l + e^). 
This corresponds to infinite 6, but in our approach we arrive at the equation 
Vg{F) = F for finite 6, which for d = 1 becomes 

9/2+x 



-J F{1~ x) dZ 



Since we know that the equation has a unique solution, it suffices to verify 
that 

Fix) = ^ + ^ 

solves it, where q satisfies 

^ ^ -21ogg 

l + q 

Via f|T3l) it can then be verified that the limit cost for finite 6 is given by 

from which the 7r^/12-result is obtained by putting g = 0. The limit cost 
(I14p for the minimum density 1 — g matching also follows in a completely 
different way from the results of [IQ]. A more streamlined derivation along 
the present lines is given in work in preparation jointly with G. Parisi. 
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3 The traveling salesman problem 



3.1 Analog of Theorem [2ll] for the TSP 

In Sections I3.2H3.7I we establish the analog of Theorem 12.11 for the traveling 
salesman problem. It follows from a theorem of A. Frieze |TT] that in the 
case d = 1, the length of the traveling salesman tour is asymptotically the 
same as the length of the polynomially solvable 2-factor problem. The idea 
is to "patch" the minimum 2-factor to a tour by replacing o{n) edges, and 
to show that this can be done at small increase in total length. The proof 
is similar to our proof of Proposition 12.141 but more complicated due to the 
global constraint of the TSP. By the concavity of the function X \^ X^^^, the 
theorem extends automatically to d > 1 but in fact, under minor changes. 
Frieze's proof works also for < d < 1. We do not discuss the details of 
Frieze's result here, but it means that we can obtain results for the TSP by 
studying the more tractable 2-factor problem. 

Our treatment of the 2-factor/TSP closely parallels the matching prob- 
lem. We focus on the differences, and omit the details where they are similar 
to those of the matching problem. 

Let Ln be the length of the minimum traveling salesman tour. The analog 
of Theorem 12.11 is 

Theorem 3.1. For every d > 1 there is a number /3TSpid) such that 

— ^Prspid). (15) 
n 

Since we divide the total length by n, Prspid) is the limit average length 
of an edge in the optimum solution. It was proved in [40j that /3t5p(1) ~ 
2.0415481864 can be expressed as 




where y satisfies 

We have computed Prspi'^) numerically, and found that 

PtspC^) ^ 1.285153753372032. 

This is consistent with the value 0.7251 given in fTEl [5U] . since the normal- 
izations differ by a factor ^/tt. 
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It is interesting to compare this value to the famous Beardwood-Halton- 
Hammersley constant [S] for the euchdean TSP in two dimensions. In [TU] this 
constant was estimated to 0.7120, but obtaining rigorous numerical bounds 
has proved annoyingly difficult. With the "euclidean" normalization (divid- 
ing by y?), our value for /3t5p(2) is 0.725070360909803, which is within 2% 
of its euclidean counterpart. It is worth pointing out that in the euclidean 
setting, the TSP does not seem to be equivalent to the 2-factor problem. 

3.2 Generalized Graph Exploration 

To carry out the analysis for the 2-factor problem, we generalize the game of 
Graph Exploration to a setting where each vertex has a nonnegative integer 
capacity. The capacity is a bound on the degree of the vertex in a feasible 
solution to the corresponding optimization problem. The matching problem 
is obtained by setting all capacities to 1, and in the 2-factor problem all 
capacities are equal to 2. In the generalized game, the rules follows: 

• The game starts at a specific vertex, and the players take turns choosing 
the edges of a walk (not necessarily self-avoiding). 

• A player who chooses an edge pays the length of that edge to the 
opponent. 

• Each player can use an edge at most once, and the set of edges chosen 
by a player must satisfy the capacity constraints, that is, each player 
can use each vertex at most a number of times equal to its capacity. 

• Moreover, Bob can use the starting point only one time less than its 
capacity. One way of thinking about this is to regard the game as 
starting by Bob entering the graph at the starting point through an 
edge coming from the outside, and that therefore Bob has already used 
the starting point once. 

• A player can, at any time, terminate the game by paying 6/2 to the 
opponent. 

3.3 The diluted flow problem 

The generalization of the matching problem to arbitrary capacities was called 
a flow problem in [IQ]. In the diluted flow problem there is a parameter 6, and 
a feasible solution is a set of edges such that no vertex exceeds its capacity. 
The cost of a solution is the sum of the lengths of its edges and a penalty for 
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vertices that are not used up to their capacity. The penalty for a vertex is 9/2 
times the difference between its capacity and the number of edges incident 
to it in the solution. 

In this setting, the proof of Proposition 12.21 goes through almost word by 
word, provided that the interpretation of G — f is that the capacity of v has 
been decreased by 1, and of course that M{G) is replaced by the cost of the 
diluted flow problem. 

3.4 Generalized Graph Exploration on the <i-PWIT 

We are led to study Generalized Graph Exploration on the rf-PWIT. An 
important difference compared to the capacity 1 case is that when the ca- 
pacities are greater than 1, it is possible to move upwards (towards the root) 
in the PWIT. Since the PWIT is a tree, a move upwards implies that the 
opponent has already used the edge in a downward move, so that after an 
upward move it is no longer possible to go back to that subtree. Moreover, 
since an upward move means that the player pays back what the opponent 
had paid to go downwards through the same edge, an upward move has the 
effect of canceling the opponent's downward move. In fact an upward move 
means that all moves played by one player in the subtree are canceled by 
moves of the opponent along the same edge in the opposite direction. 

This leads to an alternative formulation of the game: When the game 
reaches a vertex the player who is not making the next move has the right 
to forbid a number of move options equal to the capacity of v minus 1. In 
particular, when all capacities are equal to 2, it means that the player not 
to move can forbid one move option. With the alternative formulation, the 
moves that are actually carried out constitute a downward path in the tree. 
If the graph is not a tree, the alternative formulation becomes slightly more 
complicated: The effect of reversing the opponent's last move (in the original 
version of the game) is not quite the same as canceling it, since the players 
have used up one of their potential visits to a vertex that may be visited 
later on. 

Since the PWIT is a tree, the alternative formulation is correct. The 
game starts by Bob forbidding one of Alice's move options. Then Alice 
makes a move and from that point forbids one of Bob's move options, and so 
on. The similarity to the chess variants refusal chess and compromise chess 
[9] is obvious. Games where a player can forbid some of the opponent's 
move options have been called comply- constrain games in the literature on 
combinatorial games [221 [12] • I was aware of the work on these games but 
still found their connection to the TSP quite surprising. 

A valuation is now redefined as a function on the vertices of the 6'-cluster 
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that satisfies 



/(t;)=min(e/2,min2(/,-/(t;,))), 



where min2 denotes second-smallest. With the new definition, the set of 
valuations is still a lattice with a maximal element Jb and a minimal element 
/a- The crucial step in the proof of Theorem 13.11 is to show that again /a 
and /b are almost surely equal. 

In the new setting, Proposition 12.61 and Lemma 12.81 hold with obvi- 
ous modifications. We let /i^, be the measure associated to the process of 
{h, fsiyi)) with the new definition of fs- 

3.5 The branching of near-optimal play 

Again we let R be the set of paths in the ^-cluster where Bob plays optimally 
(relative to /s) and every action made by Alice is (5-reasonable. A difference 
is that now there are two different types of actions. When Alice is about to 
make a move. Bob forbids the move option that would be most preferable 
to Alice, and Alice chooses between the remaining ones. The other type of 
action is that when Bob is about to move, Alice has to forbid one of Bob's 
move options (if there are any in the ^-cluster). If Bob is about to move from 
f , and Vi and Vj are the best and second-best move options relative to /b, then 
allowing Bob to play to Vi is (5-reasonable if Ij — fB^Vj) < k — fB{vi) + S. Since 
Bob plays optimally we do not have to distinguish between Alice forbidding 
different other move options. 

We define R{k) as before, but we modify the definition of H{k) by intro- 
ducing another parameter A, assuming < A < 1. We now let 



H{k) = i^{veR{k): f^(v)<9/2} + \-i^{veR{k): fB{v) = 9/2}. (16) 



When Alice is moving, the best move option has been forbidden by Bob, 
but the analysis of Section 12.61 still goes through, and we arrive at 



We now turn to the situation when Bob is about to move. The only type 
of mistake is now if Alice forbids a move option other than the best one, 
thereby allowing Bob to play a "super-optimal" move. We begin with the 
case that Bob moves from a vertex v with JbIv) = 0/2. Then there is no 
optimal (that is, second-best) move, and the number of super-optimal moves 
is either zero or one. To upper-bound the branching we may condition on 



EH{k + l) 
EH{k) 



< max 1 - (1 - A) exp(-0"') + o(l). 




< l-(l-A)exp(-e'^)+o(l). (17) 
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exactly one point above the diagonal I — f < 9/2 in the /-/-square. The 
probability that neglecting to forbid the move corresponding to this point is 
a (5-reasonable decision is 

piv{9/2-6<l-f<e/2) 

f^vil - / < 0/2) - eM-0') 

Finally we consider the case that Bob moves from a vertex v such that 
fsiv) < 0/2. This means that there are at least two points of the process 
{h, fBivi)) above the diagonal in the /-/-square. We allow for an optimal 
(that is, second-best) move and a super-optimal move to a vertex Vi with 
fsivi) = 9/2. On the other hand we have to bound the probability of a 
super-optimal move to a vertex Vi with fsivi) < 9/2. 

Since /i^ is continuous except on the line f — 9/2, we can find an x < 9/2 
such that 

f^vix <f< 9/2) 

is as small as we please. 

If fsiv) < —X, then the probability of a super-optimal move to a vertex 
Vi with fsivi) < 6*72 is at most 

l^v(l-f< -xtf < 9/2) ^ ii^{x <f< 9/2kl < 9/2 -x) 
l^v{l - f <-x) ~ Hv{x < f kl <9/2- x) 

yjyjx <f< 9/2) n,{x < / < 9/2) 
- ii.{f^9/2) - eM-O") ' ^ ^ 

which can be made as small as we please. 

Suppose on the other hand that fsiv) > —x. Then the probability of a 
super-optimal move to a Vi with fsiyi) < 9/2 is at most 

Miv) -5<l-f< fniv)) ^ S-d- 9''-' 

IJ-vil - f < fsiv)) ~ Hv{l - f <-x)' 

By first choosing x and then choosing S we can make this too as small 
as we please. Therefore we can summarize the bounds on the expected 
contributions to H{k + 1) from the various move situations in the following 
table: 



Player to iiioa'c 


\'ortox 


Coutributioii to EH{k + 1) 


Alice moves 


iBiv) < 9/2 


l-(l-A)exp(-^'^)+o(l) 




fsiv) = 9/2 


o(l) 


Bob moves 


fB{v) < 9/2 


1 + A + o(l) 




fsiv) = 9/2 


0(1) 
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If k is even, then Alice moves from vertices at level k, and we obtain a 
recursive bound on EH{k + 1) by 

^^(^ + < max f 1 - (1 - A) exp(-^^) + "^^^ 



EH{k) - V / / ■ w> ^ 

= l-(l-A)expM^)+o(l). (19) 

When Bob moves from level k + 1, we similarly obtain 

^ ^ <max l + A + o(l),^ = l + A + o(l). (20) 



EH{k + l) - V ' " A 



Multiplying f[T9|) and fl20l) we obtain a bound on the branching effect of a 
pair of moves, one by Alice and one by Bob: 

^EH(k)^^ - ^ " "^P(-^') + ^ + exp(-^^^) + 0(1) < 1, 

for small A. It follows that if first A and then 5 are chosen small enough but 
positive, then EH{k) a.s k oo, and consequently R is almost surely 
finite. 

In the same way as in Section 12.61 it follows that there is almost surely 
only one valuation, and that E {f^i'^oot) — f\{root)) ^ as A; ^ oo. 



3.6 Integral equation for the TSP 

Let / be the unique valuation. We can state an integral equation that de- 
scribes the distribution of f{root). Let 

F{x) = P{f{root) > x) = P(at most one i such that /j — f{vi) < x), 

assuming that —6/2 < x < 6/2. The process of /j such that /(fj) > k — x is 
an inhomogeneous Poisson point process of rate dl'^~^F{l — x) on the interval 
0</<^/2 + x. Since F{x) is the probability of at most one event in this 
process, F must satisfy 

F{x) = (l + /(x))e-^(^), 

where 

i'9/2+x 

I{x) =d l^-^F{l - x) dl. 

Jo 
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It is natural to define an operator We by 



{WeG){x) 



^ + d J l'^'^G{l - x) d/ j ■ exp { -d J l^-^G{l - x) dl 

If we start with the function which is identically zero and iterate We, we 
obtain the distributions of f\{root) and f^ii^oot) for successive values of k. 
Almost sure uniqueness of valuation is equivalent to the statement that the 
sequence converges pointwise, and this in turn implies that F is the unique 
fixed point of We- 



3.7 The limit average length of an edge in the mini- 
mum tour 

In complete analogy with the case of matching, the total length Ln{0) of the 
edges in the optimum diluted 2-factor satisfies 

^-y- // ix + yr-'Fix)Fiy)dxdy, (21) 

~e/2<x,y<e/2 
x+y>0 

where F is now the fixed point of We instead of Ve- The argument of Sec- 
tion 12.91 carries over without changes. 

The quantity in the right hand side of fl2Tl) is increasing in 6, and con- 
verges to a limit l3TSp{d) as 6' — > cxd. To show that l3TSp{d) is the limit 
average length of an edge in the minimum traveling salesman tour, we have 
to do two things: First we show that Prspid) is the limit average length 
of an edge in the minimum 2-factor. The argument given in Section 12.101 
goes through without essential changes. We only need to modify the proof 
of Proposition 12.141 by introducing extra edges of another color, say blue, to 
take into account the fact that a vertex may need two more edges to complete 
the 2-factor (there are also other ways to modify the proof). At the same 
time, we may remove the assumption that n is even (this is another minor 
point that can be handled in several ways). Second, we apply the theorem 
of Frieze [TT] to conclude that Prspid) is also the limit average edge-length 
in the minimum tour. This establishes Theorem 13. 1[ 
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4 The minimum edge cover 



4.1 Another optimization problem 

Finally we turn to another optimization problem belonging to the same family 
(minimizing the total length of an edge set satisfying certain constraints). 
An edge cover is a set of edges of which every vertex is incident to at least 
one. The minimum edge cover has been studied in the bipartite pseudo- 
dimension 1 model in [13], using "finite n" combinatorial methods. The 
limit cost is the area of the union of the regions y < and x < in the 
positive quadrant. This area is equal to 1^(1)^ + 21^(1) ^ 1.456, where W 
is the Lambert H^-function defined as the inverse to the function We^ . In 
particular W{1) ~ 0.567 is the solution to x = e~^, and gives the coordinates 
of the point of intersection of the two curves y = and x = e~^. 

The edge cover problem has not been studied in the physics literature, 
but it has several interesting features and the replica-cavity method produces 
a prediction of the ground state energy with the same scheme of calculations 
as for matching and TSP. In contrast to matching and TSP, the number 
of edges in the optimum edge cover is not determined by the number n of 
vertices, but depends on the problem instance. Despite this, the edge cover 
problem is technically simpler than matching and TSP, and in particular 
allows for an explicit solution (up to numerical constants) also for d = 2. 

It turns out that the edge cover problem corresponds to a two-person 
game and that using this game, the steps of the solution can be carried 
out in analogy with matching and TSP. We try to explain how to discover 
the correct definition of this game starting from the optimization problem. 
Eventually we introduce a diluted form, but for the moment we need not 
worry about the parameter 6. 

We ask under what conditions an edge e between vertices u and v belongs 
to the minimum edge cover. For large n we can assume that e is not part of 
any short cycle, so in practice we think of e as a bridge whose removal will 
disconnect the graph into two subgraphs. Let G be the subgraph containing 
V. Alice and Bob are now advocating the options of leaving out e and of 
including e respectively, focusing on G. 

Alice, who doesn't use the edge e, simply has to find an edge cover on G. 
Bob on the other hand uses e and therefore has to find a set of edges that 
covers all vertices except possibly v. In contrast to the matching problem, 
Bob's task is not equivalent to the edge cover problem on the subgraph G — v, 
since he can use edges incident to v. 

The game therefore starts by Alice explaining how she intends to satisfy 
the constraint that Bob does not have, namely how to cover v. Possibly 
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Alice uses more than one edge from v in her cover, but she only needs to 
specify one of them in order to get rid of the constraint of having to cover 
V. Suppose Alice chooses to use the edge {v,w). Then apart from covering 
V she has fulfilled the constraint to cover w. That means that Bob is now 
behind in the sense that he has a constraint that Alice no longer has, and he 
therefore has to explain how he is going to cover w. 

Bob therefore chooses an edge from w, and there is nothing that prevents 
him from choosing the edge {v, w) that Alice is using (but the initial assump- 
tion is that Bob is using the edge {u,v), and this is the reason why in the 
game corresponding to matching. Bob does not have the right to go back to 
the starting point). If Bob chooses {v,w), the game reaches an equilibrium 
where the players' remaining constraints are exactly the same, which means 
that the game is over. Since the length of the edge {v, w) has been paid back 
and forth between the two players, the outcome of the game is zero. 

If on the other hand Bob chooses another edge {w,x), then again Alice is 
behind because she has not yet covered x and Bob has. Ahce then chooses 
an edge from x, and if that edge goes to ^; or equilibrium is reached and 
the game is over. Otherwise Alice covers a new vertex which Bob must then 
cover, and so on. 

To sum up, an "infinite 6" , or "zero temperature", version of the game is 
that Alice starts at v, players take turns choosing the edges of a path, and the 
game is over as soon as one player chooses an edge to a vertex that has already 
been visited. Since the edge cover problem has feasible solutions regardless 
of the parity of the number of vertices, the zero temperature version already 
makes sense, but to simplify the analysis we can also introduce a finite 9 
version where a constraint may be violated (a player can quit the game) at 
cost 0/2. 

If the game is played on a tree (and for finite 6 and large n it will be, in 
view of the PWIT approximation), then the rules can be reformulated: The 
players take turns choosing the edges of a downward path in the tree, and 
a player has the right to quit at cost 9/2 before their own move, or at zero 
cost after their move. The option of quitting at zero cost after their own 
move corresponds to deciding to cancel the opponent's next move. It makes 
no difference in theory if a player is required to make this decision before the 
opponent's move. 

The underlying philosophy when constructing the two-person game from 
the optimization problem is to think of an original game where Alice and 
Bob simply construct their solutions individually and then display them and 
compare the total lengths. The basic assumption is that Bob uses a particular 
edge e which Alice does not use. Assuming optimal play, it clearly doesn't 
matter who reveals their solution first, or if the players take turns revealing 
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parts of their solutions. Thus by applying principles of "strategy stealing" , 
we successively transform the game into a series of other games which are 
equivalent under optimal play. Our goal, as designers of the game, is to 
decide the outcome of the game while having the players reveal as little as 
possible of their solutions. 

Proceeding in analogy with our treatment of matching and TSP, we define 
yet another valuation concept on the PWIT by 

fiv) = max(0, min(^/2, h - f{vi))). (22) 

The crucial question of replica symmetry is whether for each 6 there is (al- 
most surely) only one valuation. If so, then f{v) and f{vi) have the same 
distribution, and we can untangle the problem from an equation of "self- 
consistency" derived from fl2^ . 

4.2 Establishing replica symmetry 

In the edge cover game, again only edges of length at most 6 are relevant. We 
define /a and fs in the obvious way, and now /^(f ) is always nonnegative. 
Therefore the relevant subset of the /-/-square is the parallelogram given by 
< / < ^/2 and < / - / < 6/2. 

The crucial question is again the branching of near-optimal play given 
that, relative to J'b, Bob plays optimally and every decision by Alice is 6- 
reasonable, where 5 is a positive number that we can choose as a function of 
d and 6. 

We want to carry out the strategy that we have used before, consisting in 
defining a certain quantity H{k) that measures the branching of R at level k. 
It turns out that we can use the same definition of H{k) as for the matching 
problem, that is, according to ([7]), and that the calculations of Section [22] 
then go through essentially unchanged. The edge cover game will terminate 
whenever it reaches a vertex of value either zero or 6/2, but to establish 
that /a = fs as in Section 12.61 we actually only have to use the fact that it 
terminates at vertices of value 6*72. 

4.3 The limit cost of the minimum edge cover 

In order to find the limit cost of the minimum edge cover, we let ^ oo. 
The fact that we can interchange the limits n — * oo and ^ oo follows from 
the results of Section and in fact the edge cover problem is considerably 
easier than matching in this respect since the analog of Poposition 12.141 can 
be established using a simple greedy scheme for completing the diluted edge 
cover. 
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We can anticipate the final results of letting ^ oo by replacing (1221) by 

f{v) = max(0, min(^/2, k - fivi))), (23) 

essentially assuming replica symmetry at zero temperature. Letting F{x) = 
P{f{v) > a^), we first notice that F{x) = 1 when x < 0. For x > 0, (125]) 
translates into 

F{x) = P(no li — fivi) < x) = P(no /j with /(fj) > U — x) 

= exp df-^F{l - x) dl^ . (24) 

Since F{1 — x) = 1 whenever Z < x, it follows that 

F{x) = exp (^-x'^ -d l'^-^F{l - x) dl^ 

= exp (^-x'^ -d (t + xy-^F{t) d?j . (25) 

4.4 The case d = 1 

When d = 1, the integral equation (I25|) takes the particularly simple form 

F{x) = e-^ ■ e-^ 

where 

poo 

A= F{t)dt 



JO 

is independent of x. It follows that A = e""^, and therefore A = W{1) and 
F{x) = W{1) ■ e"^. 

In fact the zero temperature replica symmetry for d = 1 can be verified 
directly from the dynamics of the integral operator that maps F to the 
function 

X ^ exp Fit) d?j ■ e-"". 

This operator takes A ■ to ■ e~^, and the dynamics therefore only 
consists in iterating the function A e~^. It is easily verified that this 
converges to the fixed point A = W{1) regardless of initial value of A. The 
limit cost is given by 



oo 



l-P{fi + f2>l)dl, (26) 
34 



where fi and /2 are independent and taken from the hmit distribution. In 
terms of the function F{x) = W{1) ■ e~^, we have 

P{fl + /2 > 

= (1 - W{1)) ■ W{1) • e"' + [ W{1) ■ e"^ ■ W{1) ■ e'"^ dx + W{1) ■ e"' 

Jo 

= (1 - 1^(1)) ■ W{1) ■ e"' + W{lf ■l-e-^ + W{1) ■ e-^ 

= 2W{1) ■ e"' - W{iy ■ e"' + W{iy ■ le~K (27) 

It follows that ([2nD is equal to + Vr( 1)72 ^ 0.72797 in agreement with 
the result of \X3\ for the bipartite graph. 

4.5 The case d = 2 

Also for d = 2 the problem allows a more or less explicit solution. For d = 2, 
the integral operator is given by 

exp (^-x'^ ^ ^ / x)F{t) dt^ = exp (-x^ - 2Ax - 2B) , 

/•oo 

A= F{t)dt 
Jo 



where 



'0 

and 



B= tF{t)dt. 
Jo 

This operator acts on the 2-dimensional space of A and B by 

'A\ f ■ e-2^+^'(erf (A) - 1) 

B)^ ^6-2^ + iAv^e-2^+^'(erf(yl) - 1) 

At the fixed point we must have 

B = ^e-2S - A', 

which makes it easy to numerically find A ^ 0.41079 and B ^ 0.18005. We 

get 

P{fi + /2 > = (1 - ^(0))i^(0 + [ {2x + 2A)F{x)F{l - x) dx, 

Jo 

and the limit cost is 



/■oo 

/ f ■ P(/i + /2 > /) d/ ^ 0.55872. 
Jo 
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5 Concluding remarks 



There are several slightly different random models and optimization prob- 
lems that yield to similar analysis. If the underlying graph is the complete 
bipartite graph Kn^n instead of the complete graph K^, then a slight change 
is needed in the estimates in Section 12. 4^ but otherwise the only differences 
are trivial scaling factors. For finite 6 it is possible to take the generalization 
further to various forms of near-regular graphs, but then it is not clear to 
what extent the 6 ^ oo limit corresponds to a perfect matching or tour. 

If the edge-lengths are independent and taken from some distribution 
satisfying P{1 < r) r'^ for small r, then one can show using fairly standard 
arguments that our main theorems about convergence in probability of the 
average length of the edges in the solution still hold, although the expected 
total length of the solution may not even exist. 

It is also possible to generalize our results to other optimization prob- 
lems, and this seems to lead to a number of open-ended questions. Clearly 
the results presented in this paper generalize to requiring each vertex to have 
a fixed degree other than 1 or 2, but there is also a large family of other 
problems, including minimum edge cover, where the d = 1 case can be ana- 
lyzed with the methods of [IQ]. It would be interesting to see to what extent 
one can determine their asymptotics for d ^ 1. 

There are several questions that we have left unanswered. Apart from the 
obvious question whether our results hold also for < d < 1, we have already 
mentioned the question whether the operators Voo and Woo have unique fixed 
points that are the limits of the fixed points of Vg and Wq. 

Another issue that we have not discussed is the Asymptotic Essential 
Uniqueness (AEU) of the optimum solutions. The AEU property [2] is es- 
sential for questions about the distribution of edge-lengths and the nearest 
neighbor statistics for perfect matching and the TSP. The results of [28] can 
probably be verified in the 6 ^ oo limit with our methods, but concluding 
that they are valid for perfect matching and the TSP requires interchanging 
the order of the limits n oo and 6 oo, which would be justified if the 
AEU property was established. 

We also have not discussed any algorithmic aspects. It has been recog- 
nized for a long time in the physics community that for several optimization 
problems, replica symmetry suggests that distributed iterative schemes like 
Belief Propagation and Simulated Annealing are likely to be successful. We 
have already mentioned the paper [31] which analyzes a Belief Propagation 
algorithm for minimum matching, and higher capacity problems are stud- 
ied in [1]. A natural suggestion in view of our results and in analogy with 
simulated annealing is to solve the diluted matching problem using Belief 
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Propagation and successively increase the value of 9. 
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